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Abstract 

Kinematics of geodesic flows on specific, two dimensional, curved surfaces (the sphere, hyperbolic 
space and the torus) are investigated by explicitly solving the evolution (Raychaudhuri) equations 
for the expansion, shear and rotation, for a variety of initial conditions. For flows on the sphere 
and on hyperbolic space, we show the existence of singular (within a finite value of the time 
parameter) as well as non-singular solutions. We illustrate our results through a phase diagram 
which demonstrates under which initial conditions (or combinations thereof) we end up with a 
singularity in the congruence and when, if at all, we can obtain non-singular solutions for the 
kinematic variables. Our analysis portrays the differences which arise due to positive or negative 
curvature and also explores the role of rotation in controlling singular behaviour. Subsequently, 
we move on to geodesic flows on two dimensional spaces with varying curvature. As an example, 
we discuss flows on a torus. Characteristic oscillatory features, dependent on the ratio of the two 
radii of the torus, emerge in the solutions for the expansion, shear and rotation. Singular (within 
a finite time) and non-singular behaviour of the solutions are also discussed. Finally, we conclude 
with a generalisation to three dimensional spaces of constant curvature, a summary of some of the 
generic features obtained and a comparison of our results with those for flows in flat space. 

Keywords: Kinematics, deformable media, Raychaudhuri equations, curved spacetimc and geodesic 
congruence. 

PACS numbers: 83.10.Bb, 04.40.-b, 4.20.Jb 
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I. INTRODUCTION 



Geodesic flows in curved spaces have been a topic of interest among mathematicians for 
many years. However, as far as we are aware, a detailed study of the kinematics of such 
flows, does not exist in the hterature. The kinematics of flows can either be hnked with the 
evolution of particle distribution on a manifold, or with the evolution of deformations of a 
medium itself. This aspect may be understood as follows. 



FIG. 1: Particle distributions at two time instants ti and t2 evolving under (geodesic) motion on 
a spherical surface. 

Imagine a spherical surface on which we have particles sprinkled randomly as shown in 
Fig. 1. Given some initial motion, these particles will move along trajectories governed 
by the equations of motion (geodesies, say, as a special case). In an active viewpoint, 
we might look at a region on the sphere containing a certain number of particles at 
an initial time and ask- do the particles accrete, spread out or redistribute themselves 
in an area-preserving way, as time evolves? The same question may also be posed for 
understanding deformations of a medium as well, where the particles now imply the 
points of the medium. Such investigations in a real four dimensional scenario, are of 
fundamental importance in studying accretion near black hole horizons, structure forma- 
tion and evolution of the universe at different scales, and motion of ffuid ffows to name a few. 

In [1], we restricted ourselves to flows in flat space where the evolution of each point on 
the medium was along a curve which solved the Newton's second law with damping and 
elastic forces added. What happens if we have flows on surfaces of nontrivial curvature? To 
this end, here, we look at geodesic flows on a sphere, a saddle surface as well as a surface 




3 



of varying curvature- the torus, and figure out the kinematics of flows on them subject to 
specific initial conditions. Curvature creates non-trivial effects on the flows, which need 
to be quantified. We are fortunate that in some cases, the set of equations governing the 
behaviour of flows can indeed be solved analytically. Therefore, we are able to see explicitly 
how the kinematical quantities evolve. Where analytical solutions are impossible, we resort 
to numerical methods in order to arrive at our conclusions. 

The kinematics of flows on deformable media is characterised by the kinematical 
quantities- the isotropic expansion, shear and rotation or vorticity [2] (henceforth referred 
as ESR). The deformation of the medium (or a local region in the medium) is related to the 
behaviour, with time, of the congruence of geodesies and can be characterised in terms of the 
time evolution of some deformation vector [1, 3]. The analysis of this evolution ultimately 
leads to a local, first order, nonlinear and coupled set of equations for the ESR variables, 
known in the literature, as the Raychaudhuri equations [4] . 

The Raychaudhuri equations are well-known in the study of spacetime singularities in 
gravitation and cosmology [5] -[9]. An important consequence which emerges from the anal- 
ysis of the equation for the expansion is the focusing of geodesies [10]. Moreover, the use of 
Raychaudhuri equations in the proofs of singularity theorems [11, 12] (through the notion 
of focusing of geodesies) is, today, a well-established fact [6, 7] . 

An important assumption about our forthcoming analysis may be noted here. In [1], we 
assumed that the medium (more precisely, the particles that make the medium) does not 
develop curvature due to the deformation (i.e the fiat space remains fiat). Here too, we 
assume that the deformations do not change the curvature of the medium. 

Our plan in this paper is as follows. We first derive the Raychaudhuri equations in two 
dimensional curved spaces of constant curvature. We then recall the exact solutions of 
the geodesic equations in the spherical and hyperbolic geometry, which are required as in- 
puts while solving the Raychaudhuri equations. The appearance of a finite time singularity 
depending on the initial conditions on the expansion, shear and rotation variables is inves- 
tigated analytically by including the effects of stiffness and viscous damping. The generic 
features of the solutions for different values of the parameters appearing in the Raychaudhuri 
equations and the initial conditions are clearly brought out and summarized. Some specific 
cases are also solved numerically and discussed in detail. The equations presented in this 
two dimensional case can easily be extended to curved higher dimensional spacetimes with 
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space- like sections of constant curvature (maximally symmetric spacetimes). This aspect 
of our work here is one of our motivating factors because, in future, we intend to carry 
out similar analysis in spacetimes which represent specific gravitational fields and arise as 
solutions in General Relativity. 

Following the analysis of flows on constant curvature spaces, we move on to spaces of 
varying curvature. We discuss the specific example of flows on a torus for various classes of 
geodesies on it and comment on particular features in the solutions for the expansion, shear 
and rotation. 

Finally, in the last section, we offer our comments and conclusions as well as our views 
on possible future work. 



II. GEODESIC FLOWS ON TWO DIMENSIONAL SPACES WITH CONSTANT 
CURVATURE 

A. Expansion, rotation and shear 

The deformations in a two dimensional deformable medium can be characterised in terms 
of the time evolution of the deformation vector (where i — 1,2) which connects two 
infinitesimally separated points of the medium at any time instant [1, 3]. The points at a 
later time instant move an infinitesimal distance along the respective instantaneous velocity 
directions [1]. For the purpose, the time evolution of the deformation vector for short time 
intervals can be given as, 

^ = s^(^)e + 0(A^^), (2.1) 

where the time dependent second rank tensor B'^j{t) governs the dynamics of the deformable 
medium. Further, using the equation (2.1), the second order derivative of the deformation 
vector with respect to time can be calculated as follows, 

+ B\B)j e = e- (2.2) 

In order to describe the congruence's behaviour, the evolution tensor B^j in equations (2.1) 
and (2.2) can be decomposed as the linear combination of expansion, shear and rotation in 
the following form, 

S^ = ^e5^ + <7^+a;^, (2.3) 
5 



where 9, cr'^ and a; represent, respectively, the expansion, shear and rotation. The kine- 
matics of deformations can now quantified in terms of expansion, shear and rotation (ESR) 
of the medium. The equation (2.3) can be written in the following 2x2 matrix form, 

where (j+ and Ux represent the shear parameters while ui is the rotation parameter which 
along with the expansion scalar 9, characterise the deformations of the medium [1, 3]. 

B. The evolution equations 

The evolution equations for the congruence of geodesies consist of the Raychaudhuri 
equations for the ESR variables and the geodesic equations (corresponding to the metric of 
interest) which are presented below. 

1. Raychaudhuri equations 

In order to have the evolution (Raychaudhuri) equations for ESR variable, we need to 
consider a general expression for is as follows, 

e - -R\^^v}u^i^ - K) e -PC, (2.5) 

where K^j and (3 are the stiffness and viscous damping parameters in the medium, respec- 
tively. It may be noted that the first term on the right hand side of equation (2.5) is clearly 
because of the curvature contribution and is absent while dealing with the case of flat space 
[1]. The K^j in equation (2.5) has the form given as, 

(2.6) 

kx k — k+ I 

Using equations (2.2) and (2.5) lead to, 

B) + B\B) + K) + /3B) = -R'lj^u'u"'. (2.7) 

In order to derive the evolution equations for ESR variables, let us consider a two dimensional 
metric in the following form, 

ds' = -^+r'd^', (2.8) 
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where k, — 0, +1, —1 denotes the curvature corresponding to the flat, spherical and hyperbolic 
geometries, respectively. Finally, using the decomposition of evolution tensor (2.3) in the 
equation (2.7), one obtains the Raychaudhuri equations for the ESR variables as follows. 



9 + ^9'' + 139 + 2{al + al - a;^) + « + 2A; = 0, (2.9) 

(j+ + (/3 + ^)(7+-^ + A;+ = 0, (2.10) 

dx + (/3 + ^)ax =0, (2.11) 

uj + (l3 + 9)u^0. (2.12) 



For K, = 0, these equations exactly reduce to the case of flat space time [1]. It may be recalled 
that Rij — gijR/2 in two dimensions. The analytical solutions of the equations (2.9)-(2.12) 
can be found by solving the geodesies equations corresponding to the metric (2.8). 

2. Geodesic equations 

The geodesic equations derived from the metric (2.8) are as follows, 

^ + 7T^^^' - ^ (1 - '^^') = 0' (2-13) 
(1 — Kr^j 

r(^ + 2r0 = O, (2.14) 

where the dot denotes the derivative with respect to an affinc parameter A which we shall 
identify, without any loss of generality, with the parameter t quoted earlier. The equation 
(2.14) leads to if — cjr^ where c is an integration constant. These equations can be solved 
separately for k = ±1 which are discussed below, 
(i) For K =1 : 

Let us consider r = sin(/), which transforms the equations (2.13) and (2.14) as follows, 

^ — sin0 COS0 = 0, (2-15) 

(^ + 2cot0 00 = 0. (2.16) 

The equations (2.15) and (2.16) are satisfied with = 7r/2 and ip = d (where o? is a constant) 
which represents an equatorial great circle. 
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(ii) For K =-1 : 

With r = sinh0, the equations (2.13) and (2.14) transforms in the following form, 

- sinh (f) cosh (^^ = 0, (2.17) 

(^ + 2coth(/) 0(^ = 0. (2.18) 

The equations (2.17) and (2.18) have a solution with (f> — X and Lp = f , where / is a constant. 
In order to evaluate the constants involved in the solutions of the geodesic equations, the 
normalisation condition ga/su'^u'^ = 1 leads to, 

+ r^{l- Kr^)ip'^ = {1- Kr^). (2.19) 

In view of the solution for geodesic equations (2.15) and (2.16) with (2.19), the constant 
c = d = ±1 for the case k, = 1. On the other hand c = for k = —1 for both the solutions 
considered for the geodesic equations (2.17) and (2.18). 

It is worth noting that for the corresponding case in 2 + 1 dimensions with —dt'^ added 
in the line element given in equation (2.8), one will obtain an additional geodesic equation 
corresponding to t as i — 0. Thus, t — aX + b and with the choices a — 1 and 6 = 0, we get 
back the above set of the geodesic equations. So, the generic features of geodesies, will remain 
same as they were in two dimensions. One may also note in the passing that in the case of 
cylindrical geometry with constant radius p represented by the metric ds'^ = dz^ + p^dd'^, 
the geodesic equations are obtained as z = and 'i^ = 0. In such a case, the Raychaudhuri 
equations essentially reduce to those obtained for the case of flat (Euclidean) space (see [1]). 

C. Analytical solutions of Raychaudhuri equations 

For the above mentioned choice of solutions of geodesic equations, the Raychaudhuri 
equations (2.9 - 2.12) for ESR variables reduce to the following simplified form. 



O^-e"^ + pe + 2{al + al-u'^) + K + 2k = Q, (2.20) 

&+ + {P + 0)a+-p^O, (2.21) 

&^ + {P + 9)a^-q = 0, (2.22) 

uj + {P + e)uj^O, (2.23) 
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where p — K/2 — k+ and q — —kx- The equations (2.20)-(2.23) represent the dynamics of the 
ESR variables. It may be mentioned here that for flat (Euchdean) geometry, the evolution 
equations for the ESR variables (see [1]) are recovered by putting k, — and redefining 

P = 

Using the substitution 9 = Z/Z -/3 m (2.20)-(2.23) leads to, 



-ZZ- -Z' + Z\ai + ai-u;')+i^K + 2k-^j— = 0, (2.24) 

Za+ + Za+ -pZ^O, (2.25) 

Z(Tx + Za^ -qZ = 0, (2.26) 

Zuj + Zuj = 0. (2.27) 

Solving for (T+, and u, respectively, from (2.25)-(2.27), and substituting in (2.24) followed 
by one further time differentiation yields, 

Z + 2a2Z + a^Z = D, (2.28) 

where ai = 4(p^ + g^), ^2 = K + 2k — f3'^/2, and D = —4:{pE+qF) with E and F as constants 
of integration. Now, one can write the general solution of the ESR variables as, 

. Z n pZ + E qZ + F G , , 

= ax = ^-^, u^^, (2.29) 

where G is also a constant of integration, and Z{t) is the solution of the differential equation. 
This scheme of solution of course leads to more number of constants of integration than 
conditions available from the initial values of the ESR variables. Therefore, one has to also 
ensure the satisfaction of the individual equations in (2.20)-(2.23), which yields appropriate 
number of constraint equations required for the solution of all the constants of integration 
involved. 

It is evident from (2.28) that the solution of Z{t) is dependent on the values of ai and 
a2, which are in turn dependent on the curvature, stiffness and damping parameter values. 
The solutions for all possible cases are straightforward, though may be tedious. Therefore, 
for the purpose of illustration, we consider only some special cases in detail in the rest of 
this section. The nature of solutions for all possible cases are discussed and summarized 
separately later. 
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Case 1: k — k+ — kx — and ^ = 

We start with the case without stiffness and damping. Here, we have some simplifications 
in (2.28) with a^ — l and a2 — k. The solutions for the spherical and hyperbolic geometries 
are presented below. 



Spherical geometry {k — 1) 

The exact solutions for the ESR variables in this case are obtained as follows. The 
solution for the equation for Z{t) is obtained as: 

Z{t) = (Ai + B2 + Bit) sin t + {Bi-A2- Bit) cos t - 2E. (2.30) 

From this expression, using the above-mentioned definitions for the ESR variables we obtain: 

^ ^ {A2 + Bi + Bit) cost + {B2 -Ai- Bit) sint 

{Ai + Bit) cost + {A2 + B2t)smt ' ^ ' ' 

^ 2E-{{Ai + B2 + Bit) sint + {Bi -A2- ^at) cost} 

2[(^i + Sit)cost+ (^2 + -B2t)sint] ' ^ ' ' 



F 



[{Ai + Bit) cost + {A2 + -Bst) sint] ' 



(2.33) 



C 

[{Ai + Bit) cos t + (>l2 + -B2t) sm t] ^ ' 

where, Ai, Bi, A2, -B2, E, F and G are the integration constants and will also be used 

throughout in the text for the analytical solutions of other cases as well. These can be 

defined in terms of the initial conditions on the ESR variables at t = which are represented 

as 6*0, cT+o, o"xo and ojq. The integration constants are calculated by putting back the solutions 
(2.31)-(2.34) in equation (2.9) and using the initial conditions on the ESR variables as, 

, 4 _ ^(go + 2(7+o) _ Aije^ - 2(7+o) 

'"(^o^ + Ko)-4/o' 2 ' 2 ' 

^2 = 1, £; = 0, F^AiGxo. G^AiUQ, (2.35) 

where /q = cr^ q — cJq. The constant B2 is arbitrarily chosen as unity. 



Hyperbolic geometry {k — —1) 

In this case, the exact solutions of the ESR variables can be expressed as, 

^ ^ {A2 + Bi + B2t) cosht + {Ai + B2 + Bit) sinht 

{Ai + Bit) cosh t+{A2 + B2t) sinh t ' ^ ' ^ 
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_ 2E-{{A^- B2 + B^t) sinht + {A2 - + B^t) coshQ 
'^^ ~ 2[ (Ai + Sii) cosh i + (^2 + ^2^) sinh ' ^ ' ' 



F 



[{Ai + Bit) cosht + (^2 + Bit) sinht] ' 



(2.38) 



(jr 

[{Ai + cosh t + (^2 + Bit) sinh t] ' ^ ' ^ 

The integration constants for this case remain same as given by (2.35). It may be noted 
that, in the absence of stiffness and damping, the solutions of the ESR variables in spherical 
geometry are in terms of harmonic functions, while the solutions are composed of hyperbolic 
functions in hyperbolic geometry. Therefore, while in spherical geometry a finite time sin- 
gularity is inevitable (irrespective of initial conditions), in hyperbolic geometry a solution 
singularity may be absent for specific initial conditions. This fact can be noted by looking 
at the values of t for which the denominator can go to zero in each case. For the spherical 
geometry the denominator goes to zero if: 

Ai + Bit , , 

The right hand side has value ^ at t = 0, and as t — )■ 00, we get Since the left hand 

side is not a bounded function there will always be solutions of the above equation at finite 
values of t. On the other hand, for a hyperbolic geometry, we have. 

The left hand side is now a bounded function (takes on values between —1 and +1). Thus, 
if Ai > A2 and Bi > B2, we find that there are no solutions which reach a finite time 
singularity. 



Case II: a2 — and k, — ±1 

This case involves both stiffness and damping. The choice a2 = implies 2k — 0^ /2 = ^1 
(for K = ±1) which signifies a relation between the damping parameter /3 and the parameter 
k related to the elastic force. Hence, one might be tempted to note a similarity between this 
case and that of standard critical damping (say for a harmonic oscillator). 

The solution of equation (2.28) then, consequently leads to the following ESR variables, 

^ ^ 2 7[e^*{Ei cos7t - Ai siwyt} + e~^''*{A2 sin7t - B2 cos7t}] _ 
~ e'y^ {Di sm -ft + Ci cos -ft} + e-'rt{D2 cos -ft -C2 sin -ft} ^ ' ^ 
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^ - (| + + c^'{Ai_ eos-f + Z?i siu-f} + (-^'{Ai eos-f + D2 siu-t}] 

~ 7[e^*{-0i sin 7^ + Ci cos 7^} + e-T*{L>2 cos - C2 sin 7^}] ' ^^'"^^^ 



- A;x + e'^*{Ai cos -ft + Si sin 7O + e-T*{A2 cos 7^ + B-2 sin 7i}] 
7[eT*{£)i sin + Ci cos 7^} + e-'^{D2 cos 7t - C2 sin 7t}] 



(2.44) 



G 

^ ~ 7[eT*{i:)i sin 7^ + Ci cos 7i} + e-T*{£)2 cos 7i - C2 sin 7^}] ' 
where 7 = and Ci = Ai + Bi, Di = Bi - Ai, C2 = A2 + B2, D2 = B2 - A2, Ei, E2 
and Gi are the integration constants. These constants are calculated by putting back the 
solutions in the corresponding evolution equations and using the initial conditions of the 
ESR variables as follows, 

A, = ^{2P-Q + 1), B, = i(2P + g - 1), 

A2 = ^(eo + ^Q + j), 52 = ^(-^o + 7<5 + 7), 

E = j cr+o, F = 7 CTxo, G = jujQ, (2.46) 

where P and Q are given as follows, 

P--^ [(1 + 2A;+)a+o - 2A;x<Jxo], ^ = ^ t^o - 2(^% + ^'xo - '^o)]- (2-47) 

In calculating the integration constants E and F, we have further used Z{0) — without 
any loss of generality. It may be noted that the solution of Z{t) consists of one exponentially 
growing and the other exponentially decaying harmonic mode. Thus, there is a finite time 
at which Z{t) = irrespective of the initial conditions. This happens when the denominator 
in the above equations goes to zero implying, 

tan7t= ^^"^ , !f (2.48) 

The fact that the above condition can be satisfied for some values of t confirms that the 
ESR variables can have finite time singularities. It may also be noted that a difference in 
the above-mentioned forms of the ESR variables for the k = 1 and k = —1 cases arise only 
through the relation between the allowed values of k and /3 (i.e. 2k — /3'^/2 — ^1 for k, — ±1). 
Only the value of the additive constant piece (3 changes in the functional form of 9, when 
we change k from +1 to —1. The various generic features based on the analytical solutions 
for all possible cases are presented in the next section. 
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FIG. 2: Identification of different regions in the ai-a2 parameter plane corresponding to different 
generic features of the solution of (2.28). 

D. Nature of solutions 

The generic features of the solutions of (2.28) for different values of ai and a2 are straight- 
forward to obtain from the original fourth order differential equation. To see this, substitute 
Z{t) — Y{t) + D/ai to get rid of the constant inhomogeneous term D. The equation for 
Y{t) can now be easily solved by the substitution Y{t) — e*''*. We obtain the relation 

— 2a2r'^ + ai = from which, by solving, we get: 

r = ±(^a2± ^al - ai^ . (2.49) 

Thus, the nature of the solutions crucially depends on the value of q;| — ai, the sign of a2 
and also on whether cti > or ai = 0. 

The nature of possible solutions is schematically presented in the ai-a2 parameter plane 
in Fig. 2. Based on this figure, one can comment on the nature of solutions of the ESR 
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Conditions 


Nature of Solutions 


(i) a2 > 0: Region I 

• ai > 

• cci = 

(ii) ao < 0: Region III 

• ai > 

• ai = 


Finite time singularity irrespective of initial conditions 

Initial condition dependent oscillatory solution, or finite time singularity 

Hyperbolic solutions: initial condition dependent stable solution, 
or finite time singularity 

Initial condition dependent stable solution, or finite time singularity 


a| = a\ 


Finite time singularity irrespective of initial conditions 


02 < ai: Region II 


Oscillatory solutions with one growing and one decaying mode 



TABLE I: Different possible nature of solutions of the ESR variables corresponding to different 
regions as marked in Fig. 2. 

variables, as summarised in Table I. 

It is worth noticing from the expression of a2 = k + 2k — that a large value of k 
in hyperbolic geometry can make the solutions similar to that in spherical geometry. This 
will occur for values of k such that q;| > cci > 0. On the other hand, a large value of P 
such that q;| > CKi and q;2 < 0, in a spherical geometry the flow characteristics will resemble 
those in a hyperbolic geometry. Further, with suitable values of stiffness and damping, 
the evolution of the ESR variables in a curved space (spherical or hyperbolic) can exhibit 
features of evolution in a flat (Euclidean) space without stiffness and damping. This will 
occur whenever a2 — 0, k+ — —1/2 and kx — 0. 

The nature of solutions depending on initial conditions can also be easily determined. As 
a specific interesting example, consider the case with ai — and 0:2 > 0. In this case (2.28) 
reduces to, 

Z + 2a2^ = 0, 

which has the general solution, 

Z{t) = Ai cos + Bi cos + + B2, (2.50) 

where Ai, Bi, A2 and B2 are constants of integration, which are determined from the 
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initial conditions. Now, as is evident from (2.29), the ESR variables will exhibit finite time 
singularity whenever Z — 0. This will occur whenever \A2\ < •\/2q;2(^i + Bf), which can 
now be explicitly written in terms of the initial conditions on the ESR variables. It may 
be noted that if there is no singularity, the ESR solution in this case will exhibit stable 
oscillations. 

E. Numerical examples 

In this section, we numerically investigate the nature of solutions for different initial 
conditions and parameter values (ai and in different regions as indicated in Fig. 2. 




0.5 1 1.5 2 2.5 

t 




FIG. 3: Evolution of ESR variables with parameters k = 1.0, k+ = kx = 0, P = 2.0, for (a) k = 1, 
and (b) k = —1. 

1. CKi > and a| = ai 

These parameter values lie exactly on the parabola in Fig. 2. When a2 > 0, the 
solutions of Z{t) are composed of harmonic functions with secular terms. This makes 
the solution of the ESR variables singular in finite time irrespective of the initial 



15 



conditions, as shown in Fig. 3(a). On the other hand, when 0:2 < 0, the solutions of 
Z{t) have hyperbohc functions with secular terms. In this case, it is possible to find 
initial conditions on the ESR variables for which the solutions are non-singular. One 

such case is shown in Fig. 3(b). 

2. ai — and q;| > 

In this case, the parameter values lie on the positive y-axis in Fig. 2. As discussed 
before, the ESR solutions can show finite time singularity or pure oscillatory modes, 
depending on the initial conditions. Both these cases are presented in Fig. 4(a) and 
(b) . It is interesting to note that the difference in initial conditions in the two sub-plots 
is only in the rotation initial condition (i.e., ujq). The stabilizing effects of rotation 
(also observed in [1]) is evident once again. 




(b) 



e„ = 1.0, <.^„ = 0.2, a^„ = 0.1, a>„ = 2.0 
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FIG. 4: Evolution of ESR variables with parameters k = 1, k = 1.0, 
for two different initial conditions. 



-0.5, fcx = 0, /3 = 2.0 
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(b) 



e„ = 1.0, a^„ = 0.2, a^„ = 0.1, co„ = 2.0 




FIG. 5: Evolution of ESR variables with parameters k = 1.0, k+ = 1.0, kx = 0.5, P = 2.0 for (a) 
K = 1 and (b) k. = —1, both of which lie in Region II in Fig. 2. 



3. CKi > and < oli 



These parameter values lie in Region II in Fig. 2. Hence, they will all show finite time 
singularity as discussed before. 

Two such cases have been presented in Figs. 5(a) and (b) for spherical and hyperbolic 
geometries, respectively. 

4. Oil > and q;| > ai 
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FIG. 6: Dependence of the nature of solutions of the ESR variables for two different sets of initial 
conditions with k = —1, A; = 0.5, k+ = 1.0, kx = 0.5, P = 4.0. 

We consider here parameter values in Region III in Fig. 2. The solutions of Z[t) here are 
composed of hyperbolic functions. The nature of solutions of the ESR variables are now 
dependent on the initial conditions. One can, therefore, have finite time singularity, as 
shown in Fig. 6(a), or stable solutions, as observed in Fig. 6(b). 

III. GEODESIC FLOWS ON A TORUS 

In this section, we consider kinematics of deformations on a torus which represents a 
curved medium with varying curvature. The line element for a ring torus is given as, 

ds^ = (1 + a cos (pfdtjj^ + a^dcp^, (3.1) 

where ip^cj) ^ [0, 27r] and a — ri/r2 with ri and r2 as the minor and major radius of the 
torus, respectively. 
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1. Raychaudhuri equations 



Using (2.7) and (3.1), the Raychaudhuri equations for ESR variables can now be obtained 
as follows, 



+ ^0^ + (30 + 2{al + al - u') + '""^ +2k^0, (3.2) 
2 V + X ^ a(l + acos0) ^ ^ 



a^ + {/3 + 9) <7+ - _ + fc+ = 0, (3.3) 
2a(l + a cos 0) 

ax + (/3 + 9) a^+k^ = 0, (3.4) 

d; + (;9 + ^)a; = 0. (3.5) 

One can try to solve these equations by simultaneously solving the geodesic equations cor- 
responding to the metric (3.1). 

2. Geodesic equations 

The geodesic equations corresponding to the above metric are given by, 

7 2a sin ■ ■ 
^- (l + acos0) ^'^^°' 

^^(l + acos0)sin0^,^^_ (3.7) 
a 

The first integrals of these geodesic equations are obtained as, 

^ = ^^2 ^ (3-8) 
(l + acos0)^ 



' k, (3.9) 



a^(l + a cos 0)^ 

where Zi and I2 are integration constants and their different choices lead to five different 
famihes of geodesies [13]. The normahsation gapu'^u^ — 1 for this case yields I2 — 1/a^, and 
the equation (3.7) is therefore only dependent on li. We now consider some special cases 
in order to perform a detailed analysis of the deformations as obtained by integrating the 
Raychaudhuri equations (3. 2)- (3. 5). 
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3. Kinematics of ESR for certain geodesic families 

We comment here on the nature of solutions of Raychaudhuri equations for certain 
geodesies famihes. 
Case I: = and V' 7^ 

This case corresponds to the inner equator geodesic of a torus. From equations (3.8) and 
(3.9), we have = constant and iIj = l/h- 



FIG. 7: Geodesies on a torus considered for the numerical integration of the Raychaudhuri equa- 
tions for (a) a = 0.2, (b) a = 0.5, and (c) a = 0.7. 

The Raychaudhuri equations arc given by (2.20)-(2.23) with the substitutions k = 7 = 
cos0/a(l + acos0) (which is constant in this case), p — 7/2 — k+ and q — —kx- The value 
of 7 can be positive, negative or zero. For positive (negative) value of 7, the solution can be 
obtained from the spherical (hyperbolic) solution as presented in the previous section. On 
the other hand, for 7 = (corresponding to (j) = 7r/2, 37r/2), the solution follows from the 
fiat space solutions as obtained in [1]. Therefore, the kinematics of deformations arc similar 
to the cases discussed previously. 
Case II: 7^ and ^ = 

This case corresponds to the meridians of a torus. From equations (3.8) and (3.9), we now 
have (f) — 1/a and = constant. The Raychaudhuri equations have the same structure as in 
Case I above with 7 = 7(A). It is difficult to integrate these equations analytically because 




(a) 



(b) 



(c) 
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of the variation of 7 with A. However, we can integrate them numerically. 

We now present numerical solutions of the Raychaudhuri equations corresponding to 
different values of the dimensionless parameter a. The corresponding geodesies are shown 
in Fig 7 which belong to two geodesic families on a torus [13]. Excluding the geodesies 
corresponding to the above mentioned Case I, the evolution of the ESR variables is observed 
to be oscillatory. These oscillations are basically due to periodic variations of 7 = R/2 
(where R is Ricci scalar) with 0. The frequency and amplitude of oscillations for small 
values of a is high and reduces with the increase in a. The evolution of the ESR variables is 
presented in Fig. 8 for three different values of a. Solutions may or may not exhibit a finite 
time singularity, depending on the initial values of ^-nd Uq. 

The Raychaudhuri equations for other known families of geodesies on the torus are, 
similarly, difficult to solve analytically and can be studied numerically. We do not indulge 
in discussing them further in order to avoid repeating qualitatively similar conclusions. 
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FIG. 8: Evolution of ESR variables for different values of a corresponding to the three geodesies 
shown in Fig. 7 {h = 0.3, I2 = l/a^). 
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IV. COMMENT ON EVOLUTION EQUATIONS IN THREE DIMENSIONAL ME- 
DIA 



The general expression for the evolution tensor B^j in in three dimensional media [1] is 
expressed as: 



-9 + + 



(4.1) 



^ ^0 
0^6 



le J 



Us 

-UJ2 OJi 



-OJi 





(4.2) 



where the expansion (9), shear (cr*^) and rotation (u^j) in the equation (4.1) have the following 
3x3 matrix form : 

^ /(7ii (7i2 (7i3 
+ Cri2 (^22 (^23 

^cri3 (^23 —(crii+a22) J 

It may be noted that, for simplicity, we have used the definition cr'^- = aij and uj\ — —003, 
uj\ — UJ2 and uj\ — —oui. In order to derive the evolution equations in three dimensions let 
us consider, the following three dimensional line element. 



+ r {dip + sin (pdx)- 



(4.3) 



1 — Acr^ 

This line element represents three dimensional maximally symmetric spaces of constant 
curvature which are topologically (k — I), {h — —1) or (k — O). 
The evolution equations (i.e. Raychaudhuri equations) for the expansion, shear and rotation 
can be derived following standard methods. There are nine coupled, first order, nonlinear 
equations involving 9 (expansion), an, au, ^22 ,o'i3,a23 (shear) and cui, 002-1 (^s (rotation) 
and their first derivatives. 

The geodesic equations in the above mentioned line element are as follows. 



r + 



f — r (1 — Kr ) cp — r (1 — ) sin cpx =0, 



(4.4) 



(1 - Kr^) 

r (f + 2f ip — r simp cos px^ = 0, (4.5) 

rx + 2fx + 2cot(p(px — 0. (4.6) 

For the three dimensional case, it may be noted that the equations (2.13) and (2.14) are 
modified with additional terms and wc also have a third equation (4.6). These equations 
reduce to the geodesic equations for the two dimensional case for a constant value of x- 
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We make a particular choice of geodesies for which x— constant, (j) — d, r — 1 and with 
K, — 1 (i.e. the case of spherical geometry). Further, — kSij. It is also possible to 
consistently choose Uia = (J23 = 0,u!i = LO2 = 0. With all of these choices, the evolution 
equations for the five remaining variables turn out to be, 

e+\e^ + pe + 2{al^ + al^ + al^ + aiia22 - + 3A; + 6 = 0, (4.7) 

^11 + (/3 + ^^) cTii + ^(ai\ + al^ - ul) - '^{al^ + aiia22) -2 = 0, (4.8) 

2 

(712 + (^ + + an + (722) (712 = 0, (4.9) 

(722 + (;5 + '^0) (722 + ^((721 + (jI^ - ^1) - ^(cT^, + (7n(722) - 2 = 0, (4.10) 

2 

^3 + (/3 + 3^ + (711 + CT22) t<^3 = 0. (4.11) 

These equations, as well as more general cases can indeed be solved numerically. We have 
worked through some examples (which we do not quote here). Generic features are not 
drastically different. 



V. SUMMARY AND SCOPE 



In this article, we have investigated the kinematics of flows on curved, deformable media 
by solving the evolution equations for the expansion, shear and rotation of geodesic con- 
gruences in two dimensional spaces of constant as well as varying curvature. Let us briefly 
point out the novelties that arise while considering flows on curved spaces. 

• In [1], where we looked at flows on flat space, geodesies were straight hues and the need 
for solving the geodesic equations did not arise. Here, however, this is a necessity and 

wc have looked at some simple solutions of the geodesic equations which are required 
as inputs in the evolution equations for ESR. 

• In a flat space setting, the evolution equations do not have a contribution from geomet- 
ric terms involving Ricci and Weyl. Here, though any contribution from the Weyl term 
is absent (two dimensions !), the Ricci term does indeed contribute. Thus, curvature 
effects on the nature of the flows do appear through the non-zero value of k. 
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• The analytic solutions obtained for the special cases illustrate the possibilities of the 
existence of finite time singularities in the geodesic congruence under study. We are 
able to point out, in a fairly general setting, the conditions on the parameters which 
imply the existence of finite time singularities/nonsingular solutions. The representa- 
tion of regions in parameter space where singular/non-singular solutions exist is also 
obtained and explained. 

• A major difference between flows in positive and negative curvature spaces is the fact, 
that in the latter, we do have the possibility of exclusively non-singular solutions. This 
is illustrated through the analytic solutions as well as the numerical results where we 
explicitly show the differences that arise between results in spaces of positive and 
negative curvature. 

• The stabilizing effect of rotation, as observed in the flat space analysis in [1] persists. 
We have illustrated this in our numerical solutions and plots thereof. 

• The role of stiffness, damping and the initial conditions are pointed out in each case 
with representative plots. 

• In the case of flows on the torus (varying curvature) the oscillatory nature of the 
inhomogcneous term in the Raychaudhuri equations leads to oscillations in the ESR 
variables. The nature of the oscillations (i.e. amplitude and frequency) seem to depend 
on the parameter a. Further, the solutions exhibit finite time singularities and the role 
of rotation is expected to be similar, i.e., it can inhibit the formation of finite time 
singularities. 

Finally, it is customary to ask in conclusion-what next? We have mentioned briefiy the 
equations in three dimensional spaces of constant curvature and pointed out that the results 
will be similar to those in two dimensions. It will surely be of greater value if a complete 
D-dimcnsional analysis is carried out. In two dimensions, it will be useful to look at further 
examples of flows in spaces with varying curvature- in particular we can look at higher genus 
tori and flows on them. Furthermore, a case of physical interest could be the kinematics 
of geodesic flows in black holes spacetimes, where the spaces are not usually of constant 
curvature and also involve a Lorentzian signature line element. 

We hope to report on some of the above questions in the near future. 
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